Trapping of two-component matter-wave solitons by mismatched optical lattices 
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We consider a one-dimensional model of a two-component Bose-Einstein condensate in the pres- 
ence of periodic external potentials of opposite signs, acting on the two species. The interaction 
between the species is attractive, while intra-species interactions may be attractive too [the sys- 
tem of the bright-bright (BB) type], or of opposite signs in the two components [the gap-bright 
(GB) model]. We identify the existence and stability domains for soliton complexes of the BB and 
GB types. The evolution of unstable solitons leads to the establishment of oscillatory states. The 
increase of the strength of the nonlinear attraction between the species results in symbiotic stabi- 
, lization of the complexes, despite the fact that one component is centered around a local maximum 

of the respective periodic potential. 

I. INTRODUCTION 

Optical lattices (OLs) offer a powerful and ubiquitous tool for the creation and control of various patterns in 
Bose-Einstein condensates (BECs). OLs are induced by the interference of counterpropagating coherent laser beams 
illuminating the condensate and creating an effective periodic potential for atoms A well-established fact is that 
OLs support gap solitons in BECs with repulsive interactions between atoms. Effectively one-dimensional (ID) gap 
solitons were predicted [2] and then created experimentally in BECs filling a "cigar-shaped" trap, which acted in 

I the combination with an OL induced in the axial direction. In self-attractive media, periodic potentials allow one to 
capture a soliton at a prescribed position, and also give rise to stable multi-soliton complexes. In fact, the model of a 
medium with the cubic self-attractive nonlinearity and effective periodic potential applies not only to BEC but also 
to optics, where it predicts spatial solitons in a planar waveguide with a transverse modulation of the local refractive 
index. Actually, the model was first put forward in the latter context, and fundamental solitons, as well as their 
bound states, were found in it 4]. Later, the same model was introduced in the context of the mean-field dynamics 

• in BECs Q. 

' Another topic of great interest to BECs is the study of binary condensates, which are most typically generated as 
mixtures of two different hyperfine atomic states, with opposite values of the z-projection of the atomic spin Q. The 
interactions between atoms belonging to the same and different Q hyperfine species can be controlled (including a 
possibility to reverse the sign of the interaction) by means of the Feshbach resonance. In particular, one may consider 
a binary condensate that features intra-species repulsive interactions combined with attraction between the species. 
It was predicted that the latter setting may give rise to symbiotic soliton complexes [o^], which are held together by 
the attraction overcoming the intrinsic repulsion. 

Specific types of BEC solitons have also been predicted in settings assuming a binary condensate trapped in OLs. 
In particular, the interplay of the repulsion between the two species, if combined with the dynamical effect induced 

I by the lattice (an effective negative mass of collective excitations) may give rise to ID and 2D symbiotic gap solitons 
[lOl [T]| , even in the case when the intra-species interaction is switched off. The addition of the intra-species repulsion 
helps to expand the stability region of the symbiotic gap solitons while the interplay of the OL effects, intra- 
species attraction, and repulsion between the two species gives rise to other new types of solitons complexes, such 
as semi-gap ones, with one of the components belonging to the semi-infinite gap in the OL-induced linear spectrum, 
while the other component sits in a finite bandgap JJJ . The case of the attraction between two self-repulsive species 
was recently considered in Ref. [l^, where it was demonstrated that the attraction leads to a counter-intuitive result, 
viz., splitting between gap solitons formed in each species. This effect is explained by a negative effective mass of gap 
solitons, which is one of their principal characteristic features |2|]. 

In fact, the OL potentials acting on two species in a binary BEC need not be identical. For example, if the 
actual source of the potential is not an OL, but rather a periodically nonuniform distribution of a magnetic field 
(acting in direction z and modulated along x), which couples to the atomic spin (while the basic trap is optical), 
then the effective potential will, obviously, have opposite signs for the atomic states with opposite values of the 
z-component of the spin. An issue of straightforward interest is then to consider two-component soliton complexes 
supported by pairs of mismatched lattices. In particular, effects of the mismatch on two-component ID and 2D 
gap solitons with linear coupling between the components (which corresponds to BEC loaded into two parallel 
tunnel-coupled traps) were studied in Refs. p^ . It was found that the mismatch affects the spontaneous symmetry 
breaking (SSB) bifurcation that accounts for a transition from symmetric or antisymmetric solitons to asymmetric 
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ones: the symmetric/antisymmetric states are replaced by their gMosi-symmetric/antisymnietric counterparts, but 
the bifurcation still occurs. An exception is the limit case when the phase shift between the lattices is tt (i.e., the two 
OL potentials are mutually opposite) - then, the bifurcation is replaced by a pseudo-bifurcation, in which branches of 
quasi-symmetric/antisymmetric and asymmetric solutions asymptotically approach each other, but never meet [l^. 
Effects of the mismatch between two gratings on the SSB bifurcations of optical gap solitons were also studied in the 
model of two linearly coupled parallel fiber Bragg gratings, which is an optical counterpart of the OL pair in the EEC 
model [T^ ("light bullets" in an array of Bragg gratings with alternating signs were considered in Ref. (isj). 

The objective of the present work is to study ID two-component solitons with the ordinary nonlinear attraction 
between the species, assuming that they are subject to the action of periodic potentials with opposite signs. The 
model corresponds, in particular, to the above-mentioned mixture of two atomic states with opposite values of the 
spin, that couples to a periodically modulated distribution of the magnetic field. Issues of straightforward interest are 
the existence and stability of two-component solitons. The stability problem is quite nontrivial in this setting, as, on 
the one hand, the soliton complex tends to be unstable because one of its constituent components is placed around a 
local maximum of the respective potential, but, on the other hand, the complex may be stabilized by the attraction 
between the components. We consider two varieties of the model, with either intrinsic attraction in both components 
[the respective soliton complexes are categorized as bright-bright (BB) ones], or opposite signs of the intra-species 
interaction [which corresponds to gap-bright (GB) complexes]. 

It is necessary to mention that various solutions (chiefly of the solitary wave type) in models of two-component BEG 
mixtures trapped in OLs, with attractive inter- and intra-species interactions, were considered in several works JL6, 17']. 
In most cases, the periodic potentials acting on both species were essentially identical (in-phase) jl3> although an 
example with the phase shift of n/2 between the two potentials was considered too [l3] (the latter work was chiefly 
dealing with soliton complexes of the dark- bright type). However, the settings studied in the present work were not 
considered before, to the best of our knowledge. 

The rest of the paper is organized as follows. We formulate the model in Sec. II, and summarize results for both 
types of the soliton complexes, of the BB and GB types, in Sec. III. Stability regions in a relevant parameter space 
are identified through the numerical computation of eigenvalues for small perturbations. The evolution of unstable 
solitons is explored by means of systematic direct simulations. The paper is concluded by Sec. IV. 



II. THE MODEL 



The starting point is the system of coupled ID Gross-Pitaevskii equations for the mean-field wave functions of the 
two BEG species, Ui and U2 [18,]. In the usual scaled form, the system is 

^§ - + t9M'u, + vA^)u„ (1) 

k=l 

with j,k = 1,2, where gjk = 4ajfe ^y mfl/h are effective interaction coefficients (ajk are scattering lengths of the 
respective interatomic interactions, m the atomic mass, common for both species, and the transverse trapping 
frequency), and Vj{x) are the normalized OL potentials, which are different for the two components. In accordance 
with the scaling, the numbers of atoms in the two species (their norms) are Nj = \ Uj{x)\'^ dx. As usual, we 

seek for stationary solutions to Eqs. ([T]) with chemical potentials fij as Uj{x,t) — exp(— i/ijt)Mj(x), where the real 
functions uj obey the equations 

1 ^2 2 

i^j'^j = ^2~d^ +2^9]kUkUj + Vj{x)uj. (2) 

k=l 

Localized solutions to Eqs. ^ were obtained through fixed-point (Newton-Raphson) iterations, using initial guesses 
for both components in the form of sech functions. As indicated in the introduction, the situations that we focus on 
here are of the BB and GB types, i.e., ones with self-attraction in both components (511.522 < 0), or opposite signs 
of the self-interaction (511 > and 522 < 0), respectively. 

Once the numerical solution to Eq. Q was found, its stability was investigated against infinitesimal perturbations 
{flj, bj} with growth rate A = Aioai + i'^imag, taking the perturbed solution as 

Uj{x, t) = cxp{—ifijt) [uj + aj{x) exp(At) + b*{x) exp(A*t)] , (3) 

and solving the eigenvalue problem generated by the substitution of this expression in Eqs. ([T]) and linearization. 
Sets of unstable complex eigenvalues of the resulting Hamiltonian eigenvalue problem arise in quartets, i.e., if A is 
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an eigenvalue, then so are —A, A* and —A*, hence the existence of any Areai 7^ impHes instabihty (real eigenvalues 
come in pairs, ±A, which also gives rise to instability). The eigenvalue problem was solved by using a finite-difference 
discretization [in the same way as it was used to solve the stationary version of Eq. H])]. The latter method reduces 
the problem to finding matrix eigenvalues, which can be done by means of standard numerical linear algebra packages. 

If a configuration is found to be unstable, an outcome of the instability development was concluded from direct 
simulations. To that effect, we used an integrator implemented on the basis of the fourth-order Rungc-Kutta method, 
with initial conditions taken as per exact numerically found stationary states, with the addition of weak perturbations 
seeding the instability growth. 

III. RESULTS: TWO-COMPONENT SOLITONS AND THEIR STABILITY 

A. Bright-bright sohtons 

The OL potentials acting on the two species are assumed to be of the form 

Vi{x) = Vn sin2(fca;), V2{x) ^ V02 cos^ikx), (4) 

which obviously corresponds to the phase shift of tt between the sublattices. Figures [T]|3] provide a systematic pre- 
sentation of the results, with a typical value of the wavenumber, k — tt/5. First, profiles of a pair of two-component 
solitons of the BB type, one unstable (through a pair of real eigenvalues), and the other one stable, are displayed in 
Fig. [H In either case, both components of the soliton are originally centered at a; = 0. In accordance with this, the 
component Ui, which is centered around its local potential minimum [see Eq. ^] is a stable one, while U2, whose 
original location coincides with a local maximum of the respective potential, V2{x), may be unstable. 

The evolution of the unstable soliton shown in Fig. [2] (the one pertaining to 512 — —0.04) was examined by means 
of direct simulations, as shown in Fig. [21 It is seen that the stable component {Ui, in this case), remains centered 
around the minimum of the corresponding potential well, x = 0, while the unstable component (U2) is displaced by the 
instability and ends up in an oscillatory state around a nearby minimum of its respective potential, x — n/ (2k) = 5/2. 
This is a typical scenario of the instability development observed in the present setting. 

Results of a systematic analysis of the existence and stability of the soliton family are presented in Figs. [3] and 
m Figure [3] shows a continuation of the family in 512, the cross-coupling interaction strength. At 512 — 0, the 
configuration is obviously unstable, due to the fact that one of the solitary waves (C/2) is sitting on top of a potential 
maximum (see, e.g., Ref. [Toj and references therein for a rigorous investigation of the instability). As the interaction 
between the two components gets enhanced, the stable first component creates an effective attractive potential for the 
second component, that eventually [in Fig. [3l this happens at gi2 < —0.19] overcomes the repulsive effect of the local 
potential V2, and thus makes the configuration stable. Due to this mechanism, one may speak about symbiotically 
stabilized solitary waves. As the continuation progresses towards more negative values of gi2, the second component 
eventually disappears, and the resulting single-component bright soliton is no longer affected by 1712. 

Full results are summarized in the two-parameter soliton stability diagram displayed of Fig. |4l with the grayscale 
indicating the magnitude of the unstable eigenvalue. The stability border between the linearly stable and linearly 
unstable regions is given by the curved (red in the online version) line. 

B. Solitons of the gap-bright type 

Typical examples of stable and unstable stationary states of the GB type are shown in Fig. O In this setting, the 
potentials are again taken as in Eq. ([4]). The stable component, Ui, is always of the gap type (i.e., with the repulsive 
intrinsic interaction, gn > 0), while the the self-attractive ("bright") component, U2, with 1722 < 0, may be unstable. 

The examples displayed in Fig. [5] correspond to /zi = 0.6, which actually falls in the middle of the first bandgap 
in the OL- induced linear spectrum of the first component. For this reason, the gap- type component (ui) does not 
feature conspicuous tails [10, US]- Taking a larger value of the chemical potential, for instance, ^2 = 0.88, which is 
chosen closer to the edge of the bandgap, we can generate the gap-components of the soliton with more salient tails, 
see Fig. [SI 

Simulations of the evolution of unstable complexes of the GB type, as shown in Fig. [71 reveal an essential difference 
from the instability-development pattern that was observed in the previous case (for the two-component solitons of 
the BB type). In that case, independently of the specific value of 512, the stable component (which was also Ui) 
remained centered around x = 0, while the unstable one, U2, oscillated around either of the two nearby minima of its 
respective potential, see Fig. [51 On the other hand, component U2 in the case of unstable GB solitons splits into two 
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FIG. 1: (Color online) Profiles of both components of stable and unstable solitons of the bright-bright type. The corresponding 
parameters in Eqs. ^ and Q are Vbi = V02 = 1, /Ui = —3, jj.2 = —1, 311 = —1, (?22 ~ —1, and ((?i2)at = —0.4 or 
(Si2)unst ~ —0.04, for the stable and unstable solitons, respectively. The solid and dashed lines show the profiles of the 
components and of the periodic potential, respectively. The bottom panel shows the respective spectral planes of the stability 
eigenvalues, A = Areai + iAimag, which identifies the left and right solitons as stable and unstable ones, respectively. 




FIG. 2: (Color online) The dynamics of both components of an unstable soliton complex of the bright-bright type is shown by 
means of density space-time contour plots. The parameter setting is the same as in the unstable case of Fig. [1] 



5 



1 

0.8 
Is 0.6 



to 

E 0.4 



0.2 




-0.6 -0.4 -0.2 

^12 



).8 -0.6 -0.4 -0.2 

^12 



FIG. 3: (Color online) The change of the stability of the bright-bright soliton (shown through the largest instability growth rate), 
and the variation of the norms of both of it's components, as gi2 decreases from to —1. Other parameters are Vbi = V02 ~ 1, 
Ml = -3, fi2 = -1, ffii = 522 = -1- 

FIG. 4: (Color online) The stability diagram for the soliton complexes of the bright-bright type in the plane of Vbi = V02 = Vq 
and gi2 (the stability is identified as per the maximum value of growth rate Aroai whose contour plot is shown). The other 
parameters are ni = —3, /i2 = —1, gu ~ 522 = —1. The black color means that the soliton complex is stable (Aroai = 0); the 
lighter the color, the more unstable the soliton is. The curved line indicates the stability border, as found from the numerical 
data. 




unequal parts, if the interaction between the two components is weak (see, e.g., an example for gi2 — —0.15 in the 
top two rows of Fig. [7]). The two splinters end up oscillating around different nearby potential minima, x — ±5/2. 
As | gi2| rises to intermediate values, such as 912 = —1, the unstable component U2 of the GB soliton ceases to feature 
splitting, but rather oscillates around one of the minima, while the first (stable) component, Ui, is found to oscillate 
between its original position, a; = 0, and an adjacent potential minimum, x — — 27r/(2fc) = —5 (in fact, Ui oscillates 
around the average position of component U2)- These scenarios of the instability development are displayed in the 
third and fourth rows of Fig. [7]). As the interaction between the components becomes still stronger - for instance, 
at gi2 = —1.65 - both components feature approximately synchronous oscillations (i.e., the soliton complex keeps its 
integrity, due to the strong attraction between its constituents) in the range of —2.5 < a; < 0, that is, between the 
minimum of potential Vi(2;) /maximum of V2{x) and the adjacent maximum of Vi (a;) /minimum of V2{x). 

Performing the continuation in the interaction-strength parameter, 512, we have generated characteristics of the 
GB-soliton family, a typical example of which is displayed in Fig. [51 cf. Fig. [3] for complexes of the BB type. Further, 
collecting the results for different values of V02 [while Vqi is fixed, see Eq. Q], we have produced the stability diagram 
for the GB solitons, as shown in Fig. [9l cf. its counterpart for the solitons of the BB type in Fig. |4l It is worthy to 
note that the critical value of gi2 at the stability border decreases almost linearly as Vq2 increases. 
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FIG. 5: (Color online) Profiles of both components of gap- bright soliton complexes in stable and unstable states. The parameters 
are: Voi = V02 = 1, fii = 0.6, ^12 = -0.3, gn = 1, 322 = -1, and (5i2)st = -2-0 , (912)^,-,^^ = -0.675, for the stable and 
unstable solitons, respectively. 




FIG. 6: (Color online) The same as in Fig. [5] except for fj.i = 0.88, ^2 ~ 



-0.4, and (312),, = -2.00, (912)^^,, = -0.675. 
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FIG. 7: (Color online) The evolution of an unstable soliton complex of the gap-bright type. Parameters are the same as in the 
former pictures displaying the gap-bright complexes, except that 312 = —0.15, g\2 = —1.00, and gi2 = —1.65, in the top two 
rows, third and fourth rows, and bottom two rows, respectively. 

IV. CONCLUSION 

The objective of this work was to extend the model of binary BEC mixtures to the case when the two components 
are trapped by mismatched OLs. The analysis was focused on the case of the strongest mismatch, with the sublattices 
of opposite signs. The intrinsic nonlinearity in one species was self-attractive, while in the other one it might have 
either sign; however, the inter-species interaction was always attractive. Due to the opposite signs of the two effective 
OL potentials, the bound state of two solitons could be unstable, as the second component of this state had to 
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FIG. 8: (Color online) Left panel: the largest instability growth rate for solitons of the gap-bright type. Right panel: the 
evolution of the norms of both components with the decrease of g\2 from to —2. The parameters chosen are Vbi = Vb2 = 1, 
/ii = 0.6, /i2 = -0.3, pii = 1, 522 = -1. 
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FIG. 9: (Color online) The stability diagram for solitons of the gap-bright type in the plane of Vb2 and 512. The notation is 
the same as in Fig. U Other parameters are fixed: Vbi = \, = 0.6, [12 = —0.3, gw = 1, 522 = —1. 



be centered around a local potential maximum. In both cases when the nonlinearity in the stable component (the 
one centered around the respective potential minimum) is self-attractive or self-repulsive (while the other species is 
self- attractive), we have identified stability regions for the soliton complexes, varying the strength of the inter-species 
attraction and of the OL depth. In cases when the soliton complex is unstable, direct simulations have made it possible 
to identify basic scenarios of the instability development, which amount to the spontaneous symmetry breaking (SSB) 
and resulting oscillations of one or both components. In the case of the soliton complex of the "gap-bright" (GB) 
type, the unstable ("bright") component may also potentially split into two unequal parts, which oscillate around 
different local potential minima. It is relevant to mention that it was not possible to implement a different type of the 
instability, where the gap-type component (the intrinsically self-repulsive one) would be originally set in an unstable 
position with respect to its sublattice, as a complex of such a type could never be constructed (even as a formal 
solution). For the same reason, we did not consider complexes of the gap-gap type. 

As concerns possible extensions of this work, an interesting issue is to consider the stability and dynamics of two- 
component solitons in two dimensions, that would be supported by two square lattices with the phase shift of tt 
in both directions x and y. In that case, it would be interesting to examine whether the unstable component may 
perform two-dimensional quasi-circular motion around the pinned stable one. 
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